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Abstract 

The optical field in free electron lasers can sometimes develop side- 
bands and exhibit very complex behavior. The purpose of this report 
is to develop a two mode model of the optical field in a free electron 
laser in which the optical modes are coupled by the free electrons in 
the laser. This model may be used to study the long term stability 
of the optical wave and the development of sidebands in the optical 
field. In this report the equations of motion governing the interaction 
of an electron with the optical modes will be derived. 

1. Introduction 

Initial theoretical studies of the free electron laser (FEL) used quantum me- 
chanics ([l]-[3]). However it was shown in [4] that a classical treatment of 
FELs would adequately describe most phases of operation ([5]-[10]) and re- 
sult in substantial simplification of the resulting equations describing their 
operation. 

One way to understand some of the main features of the FEL problem 
is to consider the motion of a pendulum with weak damping and a nearly 
constant tangential forcing. In a mathematical model of the FEL with a 
tapered magnetic wiggler, however, the weak damping term in the pendulum 
equation cannot be ignored because the wiggler is extremely long compared to 
its period, and one must account for the “cumulative” effects of this damping 
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over the entire length of the wiggler in order to obtain accurate solutions. 
The only way to systematically account for these higher order effects is to 
formulate a completely dimensionless problem which contains dimensionless 
parameters whose magnitudes characterize various distinguished limits. 

In a FEL frequently the optical field can develop sidebands and become 
highly complex in its behavior. The purpose of this report is to develop a 
model in which there are only two modes present in the optical field: these 
modes are coupled to one another through the electrons in the FEL. This 
model can be used to study the long term stability of the optical wave and the 
development of sidebands in the optical field. A regime of particular interest 
is that of sustained resonance in the two mode model. In this report we shall 
derive the equations of motion governing the interaction of an electron with 
the optical fields in the two mode model of a FEL. 

2.1 Equations of Motion 

We shall formulate the equations of motion as a first order system of 
ordinary differential equations. Our derivation follows that in [11]. However 
we shall be extending the derivation given there for a single optical mode 
to the two optical mode case. The convention we shall use is that upper- 
case letters will represent dimensional variables and lower-case letters will 
represent dimensionless variables. 

To determine the leading order equations, consider a single electron mov- 
ing in a helical wiggler magnetic field given by 

B„ = J5^{xcos|^Z + Ysin^z}, (1) 

where is the magnetic field amplitude and is its wavelength. The 
vectors X, Y,and Z are orthonormal and are associated with the coordinates 
X, Y, and Z, where Z is the coordinate along the direction of the magnet axis. 
We shall assume that the radiation (or signal) field is a linear superposition 
of two modes and is given by 



E, = + Ei^), (2) 

where 

E« = -E['\Z) {x sin 0,- + Y cos 0,} , i = 1 , 2 (3) 
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and 


B, = B<" + B«), 


(4) 


where 


B*'' = Z X = 1,2. 


(5) 


£<■') is 


the radiation amplitude of mode i and 






Oir 


(6) 



where is the wavelength, is the frequency, T is time, and is the 
phase of mode i, i = 1,2. 

We shall assume that the electron is relativistic, though its transverse 
velocity (the velocity orthogonal to the Z axis) is small. Let /3 be the dimen- 
sionless electron velocity (/3 is the electron velocity divided by c, the speed 
of light). We shall decompose /3 into a component parallel and a component 
perpendicular to the Z axis: 

/3 = /?||Z + /3_l, 0 < /5|| < 1. 

The equation of motion (Newton’s second law) for the electron is 

P = eEj+e/3 x(Bj-t-Bu;)-t-space charge force+radiation reaction force, (8) 

where e is the electron charge and P = 7mc/3 is the electron momentum. 7 
is the Lorentz factor, defined by 

c is the speed of light, and m is the electron rest mass. Criteria for neglecting 
the effects of the radiation reaction force and the space charge force were given 
in [12], and were derived by examining higher order terms in the equations 
of motion. These conditions will be given later in this report. 

The conservation of energy implies {pjmc^) = ec/3 • E, or 

7 = — /3 • E, (10) 

me 

To gain a rough idea of the approach to be followed, we shall assume 
initially that only a single mode in the radiation field is present. Later we 



(9) 



(7) 
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shall repeat the same calculations in more detail with two modes in the 
radiation field. 

The momentum equation ( 8 ) implies 

( 7 mc /3 ) = 7 (mc /3 ) + 7 mc /3 = 7 ^(mc /3 )(/3 • /3 ) + 7 mc /3 , ( 11 ) 

since 7 = 7 ^/? • /3 . 

Combining the momentum equation ( 11 ) and the energy equation ( 10 ), 
we conclude that 



7 mc /3 = eE, - e/3 (/3 • E^) + e/3 x (Z x EJ + e/3 x B„,. (12) 

Since /3 X (Z X Ej) = Z(/3 • Ej) — E,(/3 • Z), equation (12) may be written 
as 

7 mc /3 = e(l — /9||)Ea + eZ/3 • E^ + e/3 x Bu,. (13) 

Using the fact that 1 — ^|| <C 1 and that 

/3 -E, = /3^-E, < 1 , (14) 



we obtain the approximate equation 



/3 = (3 X Bu,. 

7 mc 



(15) 



Substituting /3 , the solution of equation (15), into the energy equation (10), 
we find 

eEs/3x . . . 

7 = sin (p, (16) 



me 



<f> = 



27TC 



1 - 



Am 1 + 

A, 272 



where </> = ( 2 ;r/A„, + 27 t/A,)Z — UjT + <f> 3 , and 

e.^3u, Au, 



2irmc^ ’ 



(17) 



(18) 



is defined as the dimensionless amplitude of the wiggler in Section ( 2 . 2 ). The 
details of the above calculations (for the case of tw’o modes) will be given in 
Section (2.3). 
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We can now define a (dimensionless) resonant electron energy 7,. for which 
^ = 0. This energy is obtained from equation (17) and has the value 

2 _ 1 + 6^ 

“ A, 2 

A free electron laser usually operates in a state for which 7 !=» 7r (see [2], [6], 
and [7]). Equation (16) indicates that electrons with (dimensionless) energy 
7r in a magnetic wiggler with wavelength A^ and field amplitude will 
spontaneously emit radiation (if jr < 0) or absorb radiation (if 7r > 0) of 
wavelength 

. _ Au,(l + fc^) 

' " 272 • 

2.2 Normalization and Nondimensionalization 





It is very important to decide on characteristic scales for the dependent 
and independent variables in the study of a system of differential equations 
by asymptotic methods. The appropriate choice of scales is usually dictated 
by the physical regime of interest represented in the system under study and 
can, in fact, be changed to suit a different regime. They key point is that 
once a set of scales is chosen all dependent and independent variables become 
dimensionless and order one in magnitude. There then is no ambiguity about 
the size or numerical magnitude of the dependent and independent variables. 
Of course, there are now dimensionless parameters present in the problem 
which are functions of the characteristic scales. It is f/ieiV orders of magnitude 
which indicate what physical regime is under study, and what is small, etc. 

We shall use Au,/27T = Lq as a characteristic unit of length and Lq/c = Tq 
as a characteristic unit of time. Let 2wfX^ = be the wavenumber of the 
wiggler. Then in terms of these units we define the dimensionless variables 

z = K,Z = ^ , t = K.cT=^. (21) 

^0 Jo 



(Note that lower case variables will always be dimensionless.) 
We now introduce the dimensionless field variables: 



bxu 



K^mc^ 



Bg , e j — 






Es, 



( 22 ) 
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where mc^ = 0.511 MeV. Using the normalizations introduced in equations 
(21) and (22), equations (8) and (10) may be written in dimensionless form 
as 



dt 



e, + /3 X (b, + bu,) + higher order terms, 



(23) 



dj 

dt 
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^ 5 - 



(24) 



In equation (23), the higher order terms include the normalized space charge 
force and the normalized radiation force. 

To characterize the order of magnitude of each term in these equations 
we need to introduce a small parameter e,0 < £ <C 1. One way to do this 
is to choose the parameter £ to be a measure of the slow variation of the 
wiggler field along the magnet axis. To do this we define e as 










1 

/ 



(25) 



where L and 1 are respectively the dimensional and dimensionless length of the 
wiggler and Lq is the characteristic length defined by the wiggler in equation 
(21). The smallness of e signifies that the magnetic wiggler produces a small 
effect which, for an electron moving along the magnet axis, lasts for a long 
time. 

Using the small parameter e we can define a “slow” variable z as 



2 = ez , 0<eCl. (26) 

The slow variable z enables us to define “slowly varying” functions of z in the 
following sense: an 0(1) function of z which varies slowly along the wiggler 
axis while maintaining its order of magnitude can be written as f(i). Note 
that df/dz = 0{e). A slowly varying function f of the dimensionless variable 
z may be written as 

/(I) = /(«) = = /(|), 

so that slowly varying means changes occur over length scales which are much 
smaller than the length of the wiggler. 

Typically the wiggler length L is several meters and the wavelength 
of the wiggler’s magnetic field is on the order of a centimeter, so e is a small 



6 



number, about 10“^. Also since the dimensionless wiggler length 1 is l/e, the 
wiggler must be considered long in the mathematical sense. This means that 
the cumulative effects of higher-order terms may be important. 

We will now denote all normalized (dimensionless) variables with an over- 
bar. The wiggler magnetic field is then written in normalized form as 

b,j, = 6u,(i){xcos z -f y sinz}. (27) 

It can be seen now from equation (19) that ^ 1 leads to 7^ ^ ie., 
extra high energy electrons are required. On the other hand, if ■C 1, 
the wiggler field is then too weak to significantly decelerate electrons. We 
conclude that the only reasonable choice is = 6u,(z) = 0(1). Because the 
operation of free electron lasers is near the resonant state defined in equation 
(19) (see references [2], [6] , [7]), and since typical values of 7r are large, 
we shall write the dimensionless energy of an electron and the dimensionless 
wavenumber of the signal in the following form consistent with equation (19): 

7 = £-"/^7, (28) 

^ = r%. (29) 

In the above equation, kg is the 0(1) dimensionless signal wavenumber 

and e" = is a second small parameter which measures the ratio Xg/Xw 
The determination of /x comes from a consideration of the physical regime 
under discussion. Recall that variables with an overbar are dimensionless 
and 0(1). 

After this preliminary analysis, let us assume again that there are two 
modes present in the radiation field. In normalized form its components 
appear as 

= {xsin^,- -f ycos 0, }, (30) 

b« = z X e», (31) 

for i= 1,2, and where 

ft = c-‘'V;\z -l) + (32) 

again for i = 1,2. Note that the total radiation field is given, as before, by 
equations (2) and (4). The parameter ^ in equation (30) is again determined 
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by the physical regime under discussion. We shall show in section [2.3] that 
the only really interesting limit for the free electron laser is ^ = 0. 

The radiation parameters and depend on the wave evolution pro- 
cess and must be determined by Maxwell’s equations. There are some open 
questions here which will be discussed later. We will see that evolves on 
a slower scale than 6^. 

There are basically three types of magnetic wigglers or undulators. There 
are those in which A^; varies, those in which varies, and those in which 
both At^ and vary. Our approach to this problem can be adapted to any 
of these wigglers. However, the work in [7] suggests that the constant A^ 
case is probably the best method for enhancing efficiency. Since it also offers 
apparent hardware advantages we will restrict our attention to this case. 

2.3 Higher Order Correction Terms 

In this section we calculate higher order correction terms for the equations 
of motion and we mention conditions derived in [12] for neglecting the space- 
charge force and the radiation reaction force. 

We begin by looking at equation governing the temporal evolution of the 
perpendicular component of the velocity, /?_l. From equations (23) and (24) 
we find 

= es + (/?||i + /fi) X (z X es + bw). (33) 

Since eg is perpendicular to 2 , we find that the last equation may be 
written 



dt 



d 



= (1 — /?||)eg + (3\\z X bw + ( higher order terms) (34) 



Now , for a single electron, ^ = /5||. We can regard this as an ordinary 
differential equation expressing z as a function of t. A little later it will be 
convenient to invert this relationship and to express t as a function of z. 

From [12], the following conditions must hold in order to neglect the 
radiation force: 



,10 



-13 



cm 



) = 



(35) 



Conditions necesssary for neglecting the space charge force are 
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(36) 



‘4 = 

Once a is determined (and it will be shown later that ^ is zero), the above 

two equations provide conditions which must hold in order to consistently 

neglect the effects of the radiation and space charge forces on an electron. In 

2 

the second equation, the ratio ^ is the ratio of the plasma frequency to the 
signal frequency. It will turn out that this ratio is a fundamental parameter 
measuring the rate of evolution of the radiation field. Consequently the 
second condition on the neglect of the radiation force plays an important 
role in the evolution of the optical field. 

Inserting the specific field quantities into equation (34) then gives 

= - /S\\b^{z) {x s'lnz - y cos z) 

— (1 — sin + t/cos0i) (37) 

-(1 - {x sin $2 y cos ^ 2 ) 

+ higher order terms, 

where b], = e’ since bs^ = z x ei' ^ 

Noting that dt = dz, we integrate equation (37) to obtain the form 

/?_L = ^ (x cos z y sin z) 

7 



(x sin z — y cos z) 


7 






3a/2+e 


^X cos 01 


- y sin0i) 


7jtp) 


3a/2+€ ^^ 3 ^ 


^X cos 02 


- y sin 02 ) 



where 0,- = e~°^'^z — i) + In equation (38), b'^ stands for the derivative 
of b^ with respect to e“z. However since b^ is a function of z = ez, it follows 
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that 



( 39 ) 



d\)^ 

d{(.°z) 






The 0(e®"/^) error in equation (38) is due to ignoring the variations of 
db^ld{t°‘^)^bi, and <f>a with z. To obtain (38) , we have also assumed that 
^40-) = O(eWJ)_ 

ie., the transverse velocity of the incoming beam can be 
ignored before it enters the magnetic wiggler. 

Now 

^ = £»/Ves = /Jj.(4''+ef') (40) 



since the component of the electron velocity parallel to the wiggler axis is 
orthogonal to the signal fields. Inserting equations (37) and (38) into the 
above expression in (40) , we find after some calculation and rearrangement 
that 



cfy 

dt 



_^2o+4 



T t(1) T t(2) 

sin(z + 0,) + sin (2 + 02 ) 

7 7 

'h _ h _ ' 

cos(z + ^,) + cos(z + 02 ) 

7 7 

_ _ 

^sin(02-^,) + f=^sin(0,-^2) 

.7^s ' 'ykl ' 



+0(c^). 

We shall now define the phases <f>i,i = 1,2 by 

<f>i = z-\-0i = z-\-c \z — ^) + 



(41) 



(42) 



We observe that since many of the parameters in the problem depend on 
z, as previously mentioned, it will be convenient to use d/dz rather than d/dt 
for the equations of motion. The velocity dz/dt can be determined from the 
relationships dz/dt = /9|| and 

;9||= (l-l/V (43) 
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Using the above equation and the expressions for 7 and /?j_ in equations 
(9) and (38) , we find after some effort, that 



^11 



Since 



it follows that 



l+c“( 








[74s" n 


(44) 


+0(e^“ 


) 






dy dz dj 

dt dt dz 


(45) 




dj 1 dj 

dz ^11 dt ' 


(46) 



We thus can combine equations (41) and (44) above to derive an expression 
for the rate of change of 7 with respect to the dimensionless distance along 
the wiggler axis, z. We find after some calculation that 



dz 



—e 






A T t(2) 

— — — Sin <pi H — — sin <p2 



l+t 



—e 






27" 

^ (fti'^COSl^l +3f^COS<^2) 







1 

TCd 

l-o 






sin 4>i H sin ^2 




7 


7 



(47) 






+0(c^+^) 






^sin(<^2-<^i) 



Equation (47) determines the evolution of 7 with respect to z, and it can 
be seen that this depends on the field quantities, and on the phases 4 >i,(t >2 
defined in (42). We now determine the evolution of these latter quantities 
with respect to z. Differentiating the expression for <f>i in equation (42) with 
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respect to z, recalling that ^ = 1//9||, and using equation (44) for 1//?||, we 
find that 



d(f>i 

dz 



= 

+0(e2“) 









'1 + 6 ^ 



7^ 



+ 



6 6 

jL I JL 

cos <pi H cos <f >2 









+ e^O)' 



(48) 



We can obtain an expression for the rate of change of (f )2 with respect to z in 
a similar fashion: 



d(j>2 

dz 




+0(e’“), 



cos 01 H T757 cos 02 









+ 



(49) 



where 0^’^' = d<f>['^/dz must be determined by using Maxwell’s equations to 
determine the evolution of the optical field. 

Note that the electron energy (7 in (47) depends on both the phases 
01,02, whose evolution is determined by equations (48) and (49). Note also 
that the evolution of these phases depends on the electron energy 7. 

If we now compare the higher order terms in the expressions for d7/dz, 
d0i/dz, and d02/dz, we see that unless ^ = 0, oscillatory solutions for 7, 0i, 
and 02 will not exist. Since the existence of oscillations in 7 implies a transfer 
of energy from the electron to the optical field, and this is essential to the 
operation of the FEL, we find that we must take ^ = 0 in equations (47), 
(48), and (49). 

At this point we have two independent small parameters: c and = p 
(as we have not yet specified a). From equation (25) c is a measure of the 
smallness of A^, compared to L, while measures the energy 7 of the 

electron (see (28)). As in [12], near the resonant state (7 m 7^), there are 
two distinguished limits corresponding to a = 0 and a = 2/3. The latter 
case is known as “slow resonance” and is substantially more difficult to solve 
than the case when a = 1, and for this reason we shall restrict our attention 
to the a = 1 case. 
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Equations (47), (48), and (49) form three coupled ordinary difFerntial 
equations for 7, <^i, and <f> 2 . The coefficients in these ODEs are slowly varying 
functions oi z = ez. Though the single mode case (see [12]) can be cast 
into Hamiltonian form, it does not seem possible to express the two mode 
equations in Hamiltonian form. One must also study how the optical wave 
interacts with and evolves according to the presence of the electron beam. 
(We have derived the equations for a single electron, but generally in the 
electron beam of a EEL there are many electrons). 

In order to maintain the coherence and monochromicity of the radiation, 
the evolution of the optical wave should be slow. Thus the leading order parts 
of the signal parameters and <j>\ must either be constant or depend on 
the prescribed slowly varying wiggler parameters. These signal parameters 
in the radiation field are determined by the wave evolution process which is 
governed by Maxwell’s equations. Though determination of the asymptotic 
behavior of the signal parameters is relatively straight forward in the one 
mode model (see [12] for example), it still is not quite clear how the wave 
evolution occurs in the two mode model. 

We leave the further details of this analysis, and the study of the three 
ODEs in (47), (48), and (49) to a future report. 
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